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S ' Abstract 

Using the expression of the fidehty for the most general Gaussian quantum states, 
the quantum fidelity is studied for the states of a harmonic oscillator interacting with an 
environment, in particular with a thermal bath. The time evolution of the considered 
system is described in the framework of the theory of open systems based on quantum 
dynamical semigroups. By taking a correlated squeezed Gaussian state as initial state, 

PL|' we calculate the quantum fidelity for both undisplaced and displaced states. The time 

evolution of the quantum fidelity is analyzed depending on the squeezing and correlation 
parameters characterizing the initial Gaussian state and on the dissipation constant and 

^ ' temperature of the thermal bath. 

1 Introduction 
> 

j1^ ■ In quantum information theory it is very important to quantify the similarity or dis- 

\^ . tinguishability of quantum states. There are several distance measures widely used to 

^D I distinguish quantum states in quantum information transfer, in particular the trace dis- 

C^^ ■ tance and quantum fidelity [Ij. Using Uhlmann's transition probability [2], Jozsa [3] 

defined the quantum fidelity of two quantum states. The fidelity quantifies how close are 
f~^ I the two states, and their separation can be measured, for example, by the Bures distance 

In the recent years there has been an increased interest in the study of Gaussian states 
of continuous variable systems used in quantum information processing [5] . Usually the 
C^ I quantum fidelity for quantum optics experiments, in particular quantum teleportation 

experiments, is calculated for pure coherent states as input states [6l[71[8]. However, in 
real experiments the input quantum states have some non-negligible degree of mixedness, 
achieved mainly due to the decoherence phenomenon which takes place during the in- 
teraction of the system with its environment. We mention here just a few of the many 
studies on quantum fidelity of Gaussian states. In Ref. [9] the authors made intense in- 
vestigations in order to understand the effects of the mixedness of states on the quantum 
fidelity achieved by experiments and the significance of the obtained results for the fidelity 
of Gaussian states. The quantum fidelity of Gaussian channels is studied in Refs. |101 [TT] 
and the fidelity for quantum information protocols is investigated in Refs. [U ^0\ [H], in 
particular for quantum cloning and continuous-variable teleportation. 

In the present paper we consider the simple system consisting of an one-dimensional 
harmonic oscillator interacting with an environment, in particular with a thermal bath. 
Our purpose is to study, in the framework of the theory of open systems based on com- 
pletely positive quantum dynamical semigroups, the influence of the environment on the 
time evolution of the quantum fidelity for Gaussian states of the considered system. The 



^ 



structure of the paper is the following. In Sec. 2 we briefly review the basic formalism for 
the calculation of the quantum fidelity for a general pair of single-mode Gaussian states. 
The time evolution of the harmonic oscillator in the theory of open quantum systems 
is described in Sec. 3. By taking a correlated squeezed Gaussian state as initial state, 
the quantum fidelity of this initial pure state and an arbitrary time mixed state of the 
considered system is calculated in Sec. 4. Conclusions are given in Sec. 5. 

2 Fidelity of Gaussian states 

For two quantum states, described by the density operators pi and p2, the fidelity F is 
defined as the Uhlmann's transition probability by [2] 

F{pi,p2)=[Tr{^P2^)^\ (1) 

The fidelity F is continuous with respect to pi and p2 and it has the following natural 
properties as a measure of the closeness of the two states: 

(i) < F{pi,p2) < 1 and F{pi,p2) = 1 if and only ii pi = p2] 

(ii) F{pi,p2) = F{p2,pi) (symmetry); 

(iii) F{pi,p2) = (V'i|P2|V'i) if one of the states is pure, i.e., when pi = \ipi){ipi\] 

(iv) F{UpiW ,Up2U^) = F{pi,p2), i.e., F{pi,p2) is invariant under a unitary transfor- 
mation U on the state space. 

Gaussian states have positive definite Wigner functions which may be interpreted 
as probability distributions in phase space. A general single-mode Gaussian state p is 
completely characterized by its first and second moments and can be represented by a 
correlated squeezed state [l2j or by a squeezed and displaced thermal state [l3]. 

The fidelity of two undisplaced thermal states was calculated in Ref. [U] and that of 
displaced thermal states in Ref. [E] . The fidelity of a large class of single- mode Gaussian 
states was obtained in Ref. p.6], using the fact that the corresponding density operators 
belong to the oscillator semigroup. 

The technical difficulty in calculating the fidelity comes from the square root of oper- 
ators. For Gaussian states this difficulty is solved, since the characteristic function of the 
square root of a Gaussian state is also Gaussian. By a successive use of the composition 
rule in the position representation. 



{x\PiP2\y) = / dz{x\pi\z){z\p2\y) (2) 

with 

{x\p\y) = exp[-(ax^ -|- dy^ + 2bxy) + lx + ky + g], (3) 

it was obtained in Ref. [16] the quantum fidelity of an undisplaced and unsqueezed 
thermal state pi and a displaced squeezed state p2. In Eq. ^ the parameters a, b, d, k, I 
and g have to fulfil the corresponding constraints imposed by the fact that the quantum 



density operator p is Hermitean, normalizable and non-negative. Denoting the amplitude 



of p2 by the row vector a 



and introducing the matrices Aj {i = 1,2) of the form 



O-pg Oipp 



(4) 



whose elements are connected with the variances Ugg, cjpp and covariance Gpq of canonical 
position q and momentum p operators through the relations a^g = 2a qq, Upp = 2app/h , 



Qpq = 2apq/h, we have [I6 



VaTs-Vs 



exp 



-Q^(Ai + A2)"^a 



where 



A = det(Ai + A2), S = (det Ai - l)(det A2 - 1) 



(5) 



(6) 



Using the property (iv) and denoting the mean amplitudes of pi and respectively p2 



by ai 



1,2), we obtain a similar formula for the quantum fidelity of two 



general Gaussian quantum states: 

2 



VaT5-Vs 



exp 



- 13'^ {A1+A2)-' 13 



(7) 



where now /3 = 02 — ai- 

If pi is a pure state, then det Ai = 1 and the fidelity ([7]) becomes 



det 



Ai + A2 



exp 



-/3T(Ai+A2)-i/?' 



(8) 



3 Master equation for the harmonic oscillator 

In the axiomatic formalism of introducing dissipation in quantum mechanics, based on 
quantum dynamical semigroups, the irreversible time evolution of the open system is 
described by the following general quantum Markovian master equation for the density 
operator p(t) [17j : 



dpjt) 
dt 






H,p{t)] + ^T.(iVjp{t),vJ] + [V„pit)vJ]). 



(9) 



Here H is the Hamiltonian operator of the system and Vj, Vj are operators on the Hilbert 
space of the Hamiltonian, which model the environment. In the case of an exactly solvable 
model for the damped harmonic oscillator, the two possible operators Vi and V2 are taken 
as linear polynomials in coordinate q and momentum p ^18. 19j and the harmonic oscillator 
Hamiltonian H is chosen of the most general quadratic form 



H = Ho + ^{qp+pq), Hq 



1 2 rncj 2 
2^P +^^- 



(10) 



With these choices the master equation ([9]) takes the following form: 



dp 

IE 



1 Z % 



R 



pp 



n' 



[?> te. p\] 



D. 



QQ I 



fl' 



-pAp.p\\ + ^{[<iAp,p]\ + \pA<i,p]])- 



[ID 



The quantum diffusion coefficients Dpp^ ^gqi ^pg ^^^ the dissipation constant A satisfy 
the following fundamental constraints: Dpp > 0, Dqg > and 



D D -D^ > -— 

^pp^qq ^pq ^ ^ ■ 



(12) 



In the particular case when the asymptotic state is a Gibbs state pq{oc) = e kr /Tre 
these coefficients have the form 1181 1191 



D 



pp 



hmuj coth ,. „ , Dn 



2kT 



' ^qq 



X — p h fiuj 

■ coth „ , ^ , Dpq = 0, 



2 muj 



2/fcT' 



(13) 



where T is the temperature of the thermal bath. In this case, the fundamental constraints 
are satisfied only if A > /x and 



(A^-/^^)coth2|^>A^ 



(14) 



The necessary and sufficient condition for translational invariance is X = p [181 [T9] . In 
the following general values X ^ p will be considered. In this way we violate translational 
invariance, but we keep the canonical thermal equilibrium state. 

From the master equation ()lip we obtain the following equations of motion for the 
expectation values ag and ap of coordinate and momentum: 



^ = -iX-p)agit) + ^apit), 
at rn 



(15) 



d(Tp{t) 



-mJ^aqit) - (A + p)(Tp{t). 



(16) 



In the underdamped case {uj > p) considered in this paper, with the notation Jl^ = uP'—p^, 
we obtain [HIE]: 



^^(0 



-\t 



(cos Ot + — sin r2i)(T„(0) H — sin Q.ta„iQ) 

il mil 



(17) 



ap{t) 



-\t 



TflUJ Li 

-— — — sin J7t(jq(0) + (cosilt — — sinr2t)(jp(0) 



(18) 



and aq{co) = (Tp(oo) 



0. 



From the master equation (jlip we can also obtain the equations of motion for the 
variances Uqq, app and covariance apq of coordinate and momentum: 



dcrqq{t) 
dt 



-2(A - p)aqq{t) + —apq{t) + 2Dqq, 



(19) 



-^ = -2(A + /i)cTpp(i) - 2moj^apq{t) + 2Dpp, (20) 

^^^ = -mco^aggit) + ^(Tpp{t) - 2Xapg{t) + 2Dpg. (21) 

Introducing the notations 



Xit) 



I' muagg (t) ^ ^ 2muDqg \^ 

app{t) 






D= I 2^ 

2D„, J 



(22) 



the solutions of these equations can be written in the form [181 US] 

X{t) = (Te^*T)(X(0) - X(oo)) + X(oo), (23) 

where the matrices T and K are given by 

/fi + in fi-rn 2lo \ /-2{x-in) o o \ 

T=— -/i-ifl^ + if] 2w,K= -2(A + iJl) (24) 

^^^^ V -w -w -2^l) V -2A/ 

and 

X(oo) = -{TK-^T)D. (25) 

Formula (j25p gives the connection between the asymptotic values (t — > cxd) of agg{t), 
(^ppit), <7pg(t) and the diffusion coefficients. These asymptotic values do not depend on 
the initial values agg{0), crpp{0), o"pg(0) and in the case of a thermal bath with coefficients 
(fT3D . they reduce to [ikfTQ] 



h hid hmuj hijj 

coth— — , crpp(oo) = ^-— coth— — , o-pg(oo) = 0. (26) 



2ma; 2kr "^^^^^ " 2 ^^^ 2kT' "'"" 

In the considered underdamped case we have 

-2Af /hi bi2 bi3\ 
Te^*T = -^^ 621 b22 623 , (27) 

V &31 &32 &33 / 

where bij {i,j = 1,2,3) are time-dependent oscillating functions given by Eqs. (3.78) in 

m- 

The relation (I12p is a necessary condition for the generalized uncertainty inequality 

ait) = agg{t)app{t) - a^^it) > — (28) 

to be fulfilled, where a{t) is the Schrodinger generalized uncertainty function. The equal- 
ity in relation (128]) is realized for a special class of pure states, called correlated coherent 
states [12] or squeezed coherent states. 



4 Calculation of quantum fidelity 

We consider a harmonic oscillator with the initial Gaussian wave function 



*(«) 



2Traqq{0) 



exp 



4^<?<?(0) 



1 - J^pM) i<l - ^,(0)? + ^^p(O)'? 



(29) 



where o"gg(0) is the initial spread, (7pg{0) the initial covariance, and o-g(O) and (Tp(0) are 
the initial averaged position and momentum of the wave packet. The initial state (j29p 
represents a correlated coherent state [12j with the variances and covariance of coordinate 
and momentum 



(7, 



m 



(0) 



2mu;' 



O", 



pp 



(0) 



hmco 



25(1 



cr. 



pq 



(0) 



hr 



2^f\ 



71' 



(30) 



where 5 is the squeezing parameter which measures the spread in the initial Gaussian 
packet and r, with \r\ < 1 is the correlation coefficient at time t = 0. The initial values 
(j30p correspond to a minimum uncertainty state, since they satisfy the equality in the 
generalized uncertainty relation ([28j): 



cr, 



m 



(O)app(O) 



a. 






(31) 



For (5=1 and r = the correlated coherent state becomes a Glauber coherent state. 

The finite temperature Schrodinger generalized uncertainty function o"(t), calculated 
in Ref. [20], has the expression 



a{t) 



fi' 



{e 



-4At 



+e 



4 
-2Atr(j + 



1 



5 + 



1 



1 



5(1 



2 coth e 



coth e + coth e 



LiJ 



/i2 cos{2VLt) 



5(1 -r2) 



5(1 -r2) 
^sin(2[7t) 2r^a;(l - cos(2(7i)) 



^ 



+ 



02^ 



02 

coth e + coth e}. 



where we have introduced the notation 



e = 



flLO 

2kf' 



(32) 



(33) 



We are interested to study the time evolution of quantum fidelity when the initial state 
is pure and the system is in interaction with its environment. If the initial wave function 
is Gaussian, then the density matrix remains Gaussian for all times (with time-dependent 
parameters which determine its amplitude and spread) and centered along the trajectory 
given by Eqs. (fTTl) and p^ representing the solutions (Tq{t) and crp{t) of the dissipative 
equations of motion [21]. We take the initial Gaussian wave function (j29p as the pure 
state 1 with the corresponding matrix Ai and the state at an arbitrary time t described 
by the density matrix p(t) as the state 2 with the corresponding matrix A2. Then we can 
apply the formula for the quantum fidelity ([8]), where the matrices Aj (i = 1,2) of the 
form (JH) are given by 



Ai 



a<?<?(0) apq{0) 
apq(O) app(O) 



aqq{t) apg{t) 

O'pqy^} 0,pp[t) 



(34) 



We obtain 



Ai + As = 2 



C^w(O) + (Tqq{t) - (o-pg(O) + (Tpg{t)) 

^ (o"pg(0) + o-pg(t)) ^ (crpp(O) + o-pp(t)) 



(35) 



and 



det ^^ ^ = ^k(O) + cT(t) + a5,(0)cTpp(t) + c7pp(0)a,,(t) - 2apg(0)cjpg(t)]. (36) 

4.1 Case of undisplaced initial state 

In the case when the initial Gaussian wave function ()29p is not displaced, one has ai = 
a2 = and the expression ([8]) simplifies: 



(37) 



det- 



Then we obtain 
F{t) = 



n 



a(0) + a{t) + agg{0)app{t) + app{0)agg{t) - 2apg{0)apg{t)] 



1 • 

2 



(38) 



Using formulas (|23p . ()26p and (|27p for calculating the variances of the coordinate and 
momentum, the expression (|38p of the quantum fidelity takes the following explicit form 
when the initial state is a squeezed state (r = 0): 



F{t) = 2{e 



-A\t 



1 



1 > 2 

6 + - ] coth e + coth e 



,-2At 



+ 



+4^il ( '^ ~ T ) sin(2r2t) coth e + 



CJ 



2-^2 



2rj2 
1 



(^ 



2 + 6^ + -. -4coth2e 



52 



^) - V(;L-coth2e) 



1 



cos(2J7t)] 



+1 + U + - coth e + coth^ e}" 2 . (39) 



In the particular case /i = 0, expression (I39p takes a simpler form: 



F(t) = 2{e; 



-4At 



1 



(5 + - I coth e + coth^ e 



g-2At ]^ 

+ ^—[2 + -5^ + 7^ - 4coth2 e 
2 0^ 



+ (2-(5^ - — j cos(2wt)] + 1+ ('^ + t) cothe + coth^e}" 



(40) 



The time evolution of the quantum fidelity (j39p as a function on the temperature, 
dissipation constant and squeezing parameter of the initial Gaussian state is represented 
in Figs. 1-3. From Fig. 1 we see that the fidelity is decreasing with time and temperature: 
monotonically, in the case of an initial coherent state (6 = 1) and non-monotonically if 
6^0. Fig. 2 represents the dependence on time and squeezing parameter. The oscillating 
behaviour of the fidelity is manifested mainly in the first part of its time evolution. Fig. 
3 represents the exponential decreasing of the quantum fidelity with time and dissipation 
constant: monotonically, in the case of an initial coherent state and non-monotonically 




15^1 




15^1 



Figure 1: Dependence of quantum fidelity F (!39l) on time t G [0, 16] and temperature T via 
C = coth j^,C e [1, 5], for u; = 1, A = 0.1, ^ = 0, r = 0. Case a) 5 = 1. Case h) 5 = 2. 

if the initial state is squeezed. The oscillating behaviour of the fidelity manifested in the 
second case is stronger for smaller values of the dissipation constant. It is interesting to 
remark that in the special case of an initial coherent state {5 = 1) and of zero temperature 
(T = 0) of the thermal bath, the quantum fidelity has a constant value in time, F(t) = 1, 
independent of the values of the parameter fi and dissipation constant A. 

The fidelity at the initial moment is evidently always F{0) = 1. For large times we get 
the following expression of the asymptotic fidelity when the initial state is a correlated 
squeezed state: 



F(oo) 



'1+U + 



6{l-r^] 



(41) 



coth e + coth e 



We notice that the asymptotic value of the quantum fidelity depends on the environment 
temperature and parameters of the initial state, but does not depend on the dissipation 
constant. If r = 0, then expression (j4ip becomes: 



F{oo) 



'1+ [6 + ^] coth e + coth^ e 



(42) 



and for a coherent state {6 = 1) we obtain: 



F(oo) 



1 + coth e 
When the temperature of the environment is T = 0, expression ()42p becomes: 

2^/6 



F(oo) 



6 + 1' 



(43) 



(44) 





Figure 2: Dependence of quantum fidelity F ( l39l) on time t G [0, 16] and squeezing coefficient 
S e [1/8,5], for u; = 1, A = 0.1, r = 0. Case a) ^ = 0, C = 1 (T = 0). Case b) /x = 0.08, 
C = 5/3 {C = coth^). 




Figure 3: Dependence of quantum fidelity F (15^ on time t G [0, 16] and dissipation constant 
A G [0, 0.3], for u; = 1, /i = 0, r = 0, C = 3 (C = cotli ^). Case a) S = 1. Case h) 5 = 2. 





Figure 4: Dependence of asymptotic quantum fidelity F{oo) (|4T1) on squeezing parameter 
6 e [1/8, 5] and temperature T via C = coth ^, C G [1, 5]. Case a) r = 0. Case b) r = 3/4. 



and respectively F{oo) = 1, if the initial state is a coherent state. For 5 ^ 1, we have 
always F(oo) < 1. 

In Fig. 4 we represent the dependence of the asymptotic value of the quantum fi- 
delity on the environment temperature and squeezing parameter of the initial state. This 
value is decreasing with increasing temperature. It is also decreasing with the squeezing 
parameter, when 5 > 1 and is increasing, when 6 < 1. 

4.2 Case of displaced initial state 

In the case of a displaced initial Gaussian wave function (p9|) , in the previous expressions 
of the asymptotic fidelity the exponential factor exp{—E) of Eq. ([8]) is present, where E 
for large times (t — > oo) has the form 



E{ 



mwo-g(O) 



oo 



1 



5(1 -r2) 



cothe -2aq{0)ap{0) 



VT 



=,+4(0) 



5 + coth e 



muj 



h 



1+ 6 + 



1 



coth e + coth^ e 



(45) 



For an initial squeezed state (r = 0) this exponential factor becomes 

1 



E{oo) 

hmuj 

and for a coherent state we obtain: 

E{oo) = 



rn^w V^ (0) I - + coth e ] + cr; (0) (5 + coth e 



1 + ( 5 + - I coth e + coth^ e 



m2u;V2(0)+(72(0) 



(46) 



(47) 
hmLo{l + cothe) 

We notice that displacing the initial Gaussian state introduces an additional decreasing 
exponential factor for the asymptotic quantum fidelity. In the particular case when the 
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temperature of the environment is T = 0, the quantum fidehty for an initial coherent 
state is given just by the exponential factor: 

^(°°) = ^"P [ 2n^ ) ■ ^''^ 

We see from Eq. (j48p that in the case of a displaced coherent state and for a zero tem- 
perature of the thermal bath, the asymptotic fidelity is strictly smaller than 1, compared 
to the case of an undisplaced coherent state, when the asymptotic quantum fidelity is 1. 
This is due to the fact that the initial mean amplitude of the density operator is decreas- 
ing exponentially to zero for asymptotically large times, as a result of dissipation, even 
in the case of zero temperature of the thermal bath (see Eqs. (fTTll . (fT8]) ). 

5 Conclusions 

Using the expression of the fidelity for the most general Gaussian quantum states, we 
have studied the quantum fidelity for the states of a harmonic oscillator interacting with a 
thermal bath. The time evolution of the considered system was described in the framework 
of the theory of open systems based on quantum dynamical semigroups. By taking a 
correlated squeezed Gaussian state as the initial state, we calculated the quantum fidelity 
for both kinds of undisplaced and displaced states. We have described the behaviour of 
the time evolution of the quantum fidelity depending on the squeezing and correlation 
parameters characterizing the initial Gaussian state and on the dissipation constant and 
temperature characterizing the environment. 

The explicit expressions of the quantum fidelity obtained in the infinite dimensional 
Hilbert space of density operators can present interest in studying the similarity of Gaus- 
sian states of quantum systems with noise and dissipation, that have a time evolution 
described by completely positive quantum dynamical semigroups. At the same time, the 
evolution of open quantum systems described by the master equation (fTTj) is equivalent 
with Gaussian quantum channels with noise and decay. For a noisy quantum channel the 
state is subject to the decoherence process, which decreases the efficacy of information 
processing, and a fundamental problem is to determine the fidelity with which the quan- 
tum information can be transmitted. The obtained results applied for a bosonic field mode 
interacting with a thermal reservoir may be therefore useful in the evaluation of the relia- 
bility of quantum information processing and communication of Gaussian states through 
bosonic Gaussian channels. Quantum cloning and continuous-variable teleportation are 
quantum processes described by bosonic Gaussian channels for which the evaluation of 
the present results for the quantum fidelity between input and output states might be 
relevant in measuring the efficiency of these processes. 
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